IFT UwB /14/2000 



THE INTEGRABLE MULTIBOSON SYSTEMS 

AND 

ORTHOGONAL POLYNOMIALS Q 



Anatol Odzijewicz*, Maciej Horowski**, Agnieszka Tereszkiewicz*** 



Institute of Theoretical Physics 
University in Bialystok 
ul. Lipowa 41, 15-424 Bialystok, Poland 



Abstract 

The strict relation between some class of multiboson hamiltonian systems and the cor- 
responding class of orthogonal polynomials is established. The correspondence is used 
effectively to integrate the systems. As an explicit example we integrate the class of 
multiboson systems corresponding to q-Hahn polynomials. 



Introduction 



The results of this paper give an effective tools to integrate some class of quantum physical 
systems. Such systems play a prominent role in quantum many-body physics, nuclear physics 
as well as quantum optics. More information about these models and their physical content 
can be found in many papers and monographs on this subject and among others one may 
consult g, Q, 0, gag] , (Kar^, p^fl , pdi ]. 

This paper is devoted to a detailed study of the mathematical structures underlying 
the multi-boson hamiltonian systems, whose dynamics is generated by the operators of the 



general form given in ( 1.1 



In the first section, it is shown, that the reduction procedure applied to the multiboson 
hamiltonian system leads to some interesting operator algebras A-jz (parametrized by a struc- 
tural function TV) which generalize in a natural way Heisenberg algebra as well as sl q (2) and 
SU q (2) quantum algebras. It was demonstrated [ Odz that the algebras An are on the one 
hand side strictly related to the integration of the quantum systems. Their description in 
terms of the coherent states leads on the other hand to the connection of these algebras with 
the theory of basic hypergeomeric series and orthogonal polynomials. 
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In Section 2 we establish the explicit relation between the spectral realization of the alge- 
bra A-ji in the space L 2 (R, da) and the representation in the space of holomorphic function. 

According to the spectral theorem, the problem of integration of the dynamical system is 
equivalent to the construction of the spectral measure da for the corresponding Hamiltonian. 
Under some assumptions on its form (see Section 3) the problem of this construction can be 
explicitly solved within the framework of q-Hahn's polynomials theory. The measure is then 
uniquely determined by the solution of ^-difference Pearson equation. 

For the sake of completeness some elementary facts from q-analysis are presented in the 
Appendix A. The Appendix B contains the proofs of some fundamental properties of q-Hahn's 
polynomials. 



1 The multi-boson systems 

This section is devoted to the detailed analysis of the symmetry properties of (N + l)-boson 
systems. The dynamics of these systems is assumed to be governed by the Hamiltonian 
operator of the form: 

H = h {a^a , . . .,a* N a N ) + g (a^a , . . . ,a* N a N ) a* . . . . , 

+a . . .a N " g Q (a a , . . . ,a N a N ) , 

where ao, ■ ■ ■ , a,N and Oq, . . . , a% are bosonic annihilation and respectively creation operators 
with standard Heisenberg commutation relations: 

[ai,a*]=5ij, [a,i,aj]=0, [a*,a*]=0. (1.2) 



The following notational convention is assumed in ( 1.1 



oj* for k t >0 

af = { 1 forh = (1.3) 

a- 1 = «r kl for h < 0. 



i k °...a k N N , k ,...,k N eZ (1.4) 



The monomial 

N ' 

can be thought of as an operator which describes the subsequent creation and annihilation 
of the clusters of the bosonic modes. 
The operator 

g (ala ,...,a* N a N ) (1.5) 

is a kind of generalization of the coupling constant. The coupling constant is replaced in our 
case by a function depending on the occupation number operators of the bosonic modes. The 
operator 

ho (ao a o, • • • ,a* N a N ) (1.6) 

can always be chosen as a free Hamiltonian being a weighted sum of the occupation number 
operators of elementary modes a^ao, • • • , a^ax'- 

h J o = uoa ao + . . . + UNa N aN- (1.7) 

The considerations at this paper will be by no means restricted to this free case however. 
The Hamiltonian under consideration (1.1) is an elementary ingredient of the most general 
Hamiltonian operator 

H= ^2 9k ...k N (ao a o, • • • ,a* N a N )ao° • . . affi , (1.8) 



2 



with the functions gk ...k N connected by the following conjugation rule. 

[9k ...k N (ao a o, ■■■ , a* N a N )]* = g-k ...-k N («S a o - h, ■ ■ ■ , a*N a N - k N ) . (1.9) 



The class of model hamiltonians corresponds to many important quantum physical 

systems. Their dynamics is generated by specific operators of the form ( |1.1| ) [ J2|, [ Karl ], 
| Kar2 ], | A-I |. For this reason the analysis of the system (|l.l| ) seems to be important and may 
shed new light on the unsolved problems of quantum physics. 

In order to analyze the quantum system described by the Hamiltonian (|l]l]), it is conve- 
nient to introduce the following operators: 

A := g (agao, • • .,a* N a N ) a^ . . .a^f. (1.10) 

and 

N 

Ai = A* := ^aijtfdj , (1.11) 

3=0 

where i = 0, 1, 2, . . . , N. One assumes that real (N + 1) x (N + l)-matrix 
a = (ctij) satisfies the conditions 

deta / (1.12) 

N 

^Oiijkj = 8 i0 . (1.13) 

The operators Aq,A\,...,An,A and A* do satisfy the following commutation relations 

[A ,A]=-A, [A ,A*]=A* (1.14) 
[A,A i ] = Q = [A\A i ] (1.15) 



for i = 1, . . . N, and 

for i,j = 0,1,...,N. 
One has in addition 



[A i ,A j ]=Q (1.16) 



A* A = \g (apao - k , . . . ,a* N a N - k N )\ 2 V ko (ao a o - k ) ■ ■ -T>k N (a* N a N - k N ) (1.17) 
-4-4* = \go{a* a ,...,a* N a N )\ 2 V ko (a* a )...V kN (a* N a N ), (1.18) 

where V ko (ao a o) > ••• ^k N (ajv ^) are polynomials: 

( a k (a*) k = (a*a + 1) . . . (a*a + k) fork>0 
V k (a* a) := a k a~ k = | 1 fork = (1.19) 

1 (a*)~ fc a- fc = a*a(a*a- 1) . . . (a*a - fc + 1) fork<0. 

The operators and ^4^4* are diagonal in the standard Fock basis 

|n , m, • • • , njv) = . . l , Kf° . . . (a* N ) nif |0> , (1.20) 

where (n , m, . . . , n N ) € : = (Z+ U {0}) X . . . X (Z+ U {0}) (JV + 1 times). 

Let us note that operators j4o, . . . , An are unbounded. Whether A and j4* are bounded 
or not depends on the choice of the structural function go. All of them are defined on the 
common domain D spanned by finite linear combinations 

\v) = ^ (Ho,ii,...,iN \ n io > n h ' • • • ; n iN ) (1-21) 
(i ,ii,...,i N )eF 
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of the Fock basis elements, where F is some finite set of multi indices. 

Identifying the 1-dimensional spaces C |no, n\, . . . , nj\r) with the elements of Z^ +1 we 
obtain the action of A and A* (and their natural powers) on Z^ +1 . It is easy to see that the 
orbits of these actions are located on one dimensional lines which are parallel to the vector 
(ko,h,...,k N )eZ N + 1 . 

If the function go of ( |1.1| ) is regular and nonvanishing in all points of Z^ +1 , we can make 
some simple but useful observations. The first one is that if ^ < for some i,j G {0,1,... , N} 
then for any element \uq, ni, . . . , njv) of the Fock basis there exists M G N such that 

A M \no, ni, . . . , njv) = and {A*) M |no, »i, • • • , n^) = 0. (1-22) 

This means that orbits of A and A* in Z^ +1 are finite. In the opposite case the orbits contain 
infinitely many points. 

The second observation is that A |no, rai, . . . , rtjv) = if and only if there exists i E 
{0, 1, . . . , N} such that ki > and n, E {0, 1, . . . , fcj — 1}. Hence each orbit of >1 and A* in 
Z^ +1 has exactly one vacuum (the point anihilated by A). It is then natural to introduce the 
parametrization of the Fock basis which is in agreement with the above orbit decomposition 
of 

Replacing the occupation number operators a^ao, . . . , a* N ai\r by the operators Aq, A%, . . . ,An 
m dLVfj) and (jOf) one obtains: 

A* A = g(A -l,A u ...,A N ) (1.23) 
AA* = G(A ,A 1: ...,A N ) (1.24) 

with the function Q uniquely determined by go, polynomials Vk , ■ ■ ■ , Vk N and the linear map 

The Hamiltonian (|1.1| ) can be reexpressed in terms of ( 1.1C| ) and (1.11) in the following 
form: 

H = H (A ,A 1 ,...,A N ) + A + A* (1.25) 

It is clear that it admits N commuting integrals of motion A\ , . . . , Ajsf. [A{ , H ] = for 
i = 1, . . . , N. They commute with the operator too. This maximal system of commuta- 
tive observables is diagonalized in the Fock basis and the eigenvalues of Aq,Ai,...,An on 
|no, ni,... , njv) are given by 

N 

\i = J2 a v n i i = 0,l,...,N. (1.26) 

3=0 

The eigenvalues (Ao, Ai, . . . , Ajv) form a discrete convex cone A + C R N+1 . It is spanned 
by the columns of the matrix (ajj) with entries from (Z+ U {0}): A + = a ^Z^ +1 ^j . The 

sequences (Ao, Ai, . . . , Xn) G A + will be used as a new parametrization {|Ao, Ai, . . . , Xn)} of 
the Fock basis elements. 

In order to integrate the system (|l.l[) one can reduce it to the eigen subspaces 7ixx...x N C H 
spanned by the eigenvectors |Ao, Ai, . . . , Ajv) with fixed Ai, . . . , Ajy. This subspace is invariant 
with the respect to the algebra A re d , which is generated by the operators A, A*, Aq . These 
operators do satisfy the following relations ( 1.23j ),( |i~24 ),( 1.14|) : 

[A ,A] = -A, [A ,A*]=A*, (1.27) 
A* A = g{A -l,X 1 ,...,X N ), (1.28) 
AA* = G(A ,X 1: ...,X N ). (1.29) 
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Hence one can conclude that the problem of integration of the system (|T]l|) amounts to the 
integration of the system described by the reduced Hamiltonian: 



H red = H (A , Ai, . . . , Aat) + A + A* (1.30) 

being an element of the algebra A Te d- 

The orthonormal basis of the Hilbert subspace 7i\ x ,..\ N is formed by the vectors | Ao, Ai, . . . , \n 
with Ao such that (Ao, Ai, . . . , \n) G A+. From (|L2^-|L2^) it follows that 

^4.0 | Ao, Ai, • • • , Ajv) = Ao |Ao, Ai, . . . , Ajv) (1-31) 
A | Ao, Ai, • • • , Ajv) = \/Q (Ao — 1, • • • , A at) I A — 1, Ai, . . . , Ajv) (1.32) 
A* | Ao, Ai, . . • , Aat) = \/ Q (Ao, Ai, . . . , Aat) |Ao + 1, Ai, . . . , Aat) . (1.33) 

Let us note here that if (Ao, Ax, ... , Xn) G A + then either (Ao — 1, ... , Aat) € A + 
((A + 1, . . • , Aat) € A+) or 

A|A ,Ai,...,Aat) = (A*\X o ,X 1 ,...,X N )=0). (1.34) 

Because of ( |L28 - L~29| ) the conditions Ql.34| ) are equivalent to 

g(A o -l,Ai,...,A iV ) = (G(X ,X 1 ,...,X N ) = 0). (1.35) 

Since A + is a discrete convex cone we can easily see that the representation of A re d in 
7~L\i...\pi splits into irreducible components. These components are generated out of vacuum 
(or antivacum) states | Ao, Ai, • • • , Aat). The vacuum states are parametrized by the solutions 
Ao of the equations ( 1.35| ). In all cases under consideration the operator Aq is diagonal while 
the operator A is weighted unilateral shift operator. One does not exclude the case when the 
irreducible representations generated by | Ao, Ai, • • • , Aat) are of finite dimension. 

We can give now the following 

Proposition 1.1 If the structural function Q is regular then: 

1) dim7ix 1 ...x N < oo if and only if there exists a pair i,j S {0, 1, . . . , N} such that ^ < 0. 
if ko > then the equation A \ Ao, Ai, . . . , Aat) = is solved by 

I 1 N 

^ = jr o ~ (i-36) 

where are matrix elements of a" 1 , I E L := |0, -kn, -/co, ■ • • , — -&o| and k is the 
biggest common divisor of the numbers (ko, fci, ■ ■ ■ ,kN). Moreover Hx 1 ...x N splits onto 
the irreducible components 

^..^ = 074..^ (1-37) 

where 7v Xl Ajv are generated by A re d out of the states \Xqj, Ai, . . . , Aat). 
3) dimHAi,...,Ajv = 00 */ an< ^ on ^ i/drmTi^ Xn = oo for all I. 

Proof : follows immediately from the observations above. Irreducibility of Tv^ x N ^ s a 
consequence of the fact that | Ao,/ , Ai, . . . , Aat) is the unique vacuum in this space. QED 

The construction presented above generalizes the one of | Karl ],| Kar2 |, where the models 
with Q being a polynomial of special type are considered. In these papers the approximative 
methods of integration of the models were used. 
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As it was mentioned in the Introduction, our aim is to study some integrable family 
of Hamiltonians ( 1.30| ). Therefore, instead of the operator Aq, we will use the Hermitian 
operator 

Q. = q A -X ^ ( 13 g) 

where < q < 1. 

The structural relations ( 1.27 - L29| ) acquire the following form in terms of the operators 
A, A* and Q 



QA* = qA*Q, qQA = AQ 
A* A = 1Z(Q) 
AA* = 7Z(qQ) , 



where structural function 1Z is given by 

logQ 



n{Q) = g 



logg 



+ -Vz — 1, Ai, . . . , Aat 



(1.39) 
(1.40) 
(1.41) 



(1.42) 



The function 1Z takes positive values in all points {q n }^ = i and moreover 1Z{1) = 0. The 
algebras of the type above were analyzed in | Odz| . The reduced Hamiltonian ( |1.30 ) can be 
rewritten as 

H red = V{Q)+A + A\ (1.43) 



with the function T> given by 

V (Q) = H 



logg 



+ A ,«, Ai, . . . , Ajv • 



(1.44) 



by the use of (|L38p. 

The analysis below is restricted to infinite-dimensional case only. The discussion of finite 
dimensional case will be presented in the separate paper. 



2 Spectral and holomorphic representations 

Let A-ji be the operator algebra generated by the operators A, A* and Q. In this section 
we describe two natural and, important from physical point of view, representations of the 
algebra An- 

The first representation, which we will call a holomorphic one, is related to the coherent 
states \z) , \z\ < TZ (0), of the anihilation operator A £ An 

A\z) = z\z). (2.1) 

Let us consider the case when the orthonormal basis of the Hilbert space 7i re d '■= \ N 

\n):=\X 0jl + n,X 1 ,...,X N ), nENU{0} (2.2) 

is infinite. The basis vectors are generated by the operator A* out of the vacuum state |0). 
The action of the algebra generators on the vectors of this basis is: 

Q\n) = q n \n) (2.3) 
A\n) = ^lZ{q n ) \n- 1) (2.4) 
A*\n) = ^/lZ(q n+1 ) \n+ 1) . (2.5) 
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The coherent state \z) is thus given by 

oo n 

\z):=Y—= Z |n), (2.6) 

where z G B = {z G C : \z\ < 1Z(0)}. The states \z), z £ D form a linearly dense subset in 
Tired- Therefore, the map 

I(v)= f (v \z) (2.7) 

where v G Ti re d an d z G B , is an antilinear and one-to-one map of the Hilbert space TC re d into 
the vector space O (B) of holomorphic functions on the disc B . It was shown in Odz| that 



the image / (J~L r ed) is isomorphic to the Hilbert space L 2 (B, dfin) containing holomorphic 
functions on B which are square integrable with respect to the measure 

log x 

dn n (z,z) = 1 . t j- ^7 ( - ) lim —, ° ,g ^°° d q xdtp, (2.8) 

2?r Exp n (x) (1 - q) {q; q)^ \x J a^oo {ax^q)^ 

where z = y^e 11 " , d q x is the Jackson measure (see Appendix A), dip is the Lebesque measure 



on the circle S . The ^-exponential function Expn (see [OdzQ and the function 7 are defined 
by 

00 x n 

Exp n (x) := (z\z) = > — — . , (2.9) 

F y ' X 1 1 ^n{q)...-R{q n y V ; 

n=0 



00 / n 



with 



T (*) := E ( E ^ **■ ^ 

n=0 \fc=0 W ' / 

n{q;q) k :=n{q)n{q 2 )...n{q k ), (2.11) 

(g;g) fe :=(l-g)(l-g 2 ). ••(!-?*), (2.12) 

00 

(x;q) 00 :=\{{l-q k x). (2.13) 

fc=0 

Let us remark here, that the existence of the morphism is equivalent to the existence of 
the resolution of the unity of the type: 

\z)(z\dpji(z,z) = l. (2.14) 

By holomorphic representation of the algebra A-ji we will understand the representation 
in the Hilbert space 1?0 (B, dfin). Straightforward calculation shows that: 

Aip(z) = dnip(z) (2.15) 

A*<p(z) = zip(z) (2.16) 

Q<p(z) = <p(qz) (2.17) 

H red tp(z) = (V{Q) + z + dn)ip{z), (2.18) 

where ip G L 2 (B, d^L-ji) an d dn is a 7£-difference operator given by 

dTVp{z) := K(qQ)d cp(z), (2.19) 

Wz) := ^~^°\ (2.20) 
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see [Odz]. In the cases with 



TZ(x) 



x 



(2.21) 



dn is a g-derivative d q , and the standard derivative jj^ is obtained from d q in the limit q — ► 1. 
The quantum algebra of 1Z given by (2.21) is a q -deformation of the Heisenberg algebra. 
Hence the analytic realization of A-r introduced above is a natural generalization of the 
Bergman-Fock-Segal representation of the Heisenberg algebra. 

The second representation of A-jz is related to the spectral measure of a selfadjoint exten- 
sion of the Hamiltonian ( [1,43 ). The action of H re d on the elements of the orthonormal basis 
{|n)}^ =0 is given in terms of three-diagonal (Jacobi) matrix 



Hred\n) = ^TZ (q n ) \n - 1) + V (q n ) \n) + (q n +±) \n + 1) 



(2.22) 



We will call this matrix the Jacobi matrix of the operator H re d- 

The operator H re d is symmetric and its domain Djj red contains all finite linear combina- 
tions of the basis elements {|^)}^L - The theory of such type operators is strictly related 
with the theory of orthogonal polynomials [ A.-G ] , [jAfl , []ChJ , [ 5u | . 

Let JC W denote the deficiency subspace of H rec [ for w 6 C and Ira lo ^ 



ICu, ■= {(Hred ~ Ul) D Hr 



The deficiency indices (n + ,n_ 



n + = dim }C U , for Imuj > 
n_ = dim JC U , for Imuj < 



(2.23) 



(2.24) 
(2.25) 



of the operator H re d are (0,0) or (1, 1). In order to show this property, one should observe 
that \v) 6 ICu if and only if 

H; ed \v)=u\v). (2.26) 
where H* d is the Hermitian conjugate of H Te d- The vector 



(w) \n) G Hred, 



(2.27) 



n=0 



solves ( [2.26 ) if and only if coefficients P n (lo) do satisfy the three term recurrence equation 

(2.28) 



uP n (w) = y/n^P^i (uj) + V (q n ) P n (u) + y /n(q n + 1 )P n+1 (lo) 
n £ N, with the initial conditions 

u-V(l) 



P (u) = l, Pi (lo) 



U(q) 



and 



J2\Pn(L0)\ 2 < +00 

n=0 



(2.29) 



(2.30) 



Hence n + and n_ are equal to or 1. 

Because every P n (lo) is a real polynomial of degree n of the complex variable lo one has 
n + = n_. 

Following |A| we will call the Jacobi matrix of H re( i to be of the type D or C if the 
deficiency indices of H red are (0,0) or (1,1) respectively. 
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Proposition 2.1 i) If 

V — = +00 (2.31) 

t/ien t/ie operator H re< i has deficiency indices (0, 0) . TTiis is equivalent to its essential 
selfadjointnes. 

ii) If the set of the coherent states \z) of the annihilation operator A is parametrized by the 
disc D of finite radius TZ(0) < +00 then H re d is essentially self adjoint. 

Hi) If the deficiency indices of H re d are (1,1) then coherent states \z) of A exist for any 
zeC. 

Proof: 



i) Let {Q n (uj)}™ =0 be an another solution of the recurrence (2.28) with the initial condi- 
tions given by : 

Qo H = 0, Qi (w) = 



1 



One then has 



Pk-i M Q k M - -Pfc M Qfe-i (w) 



where G N. Applying Schwartz inequality to (|2.33| ), one finds 



<2 £|p„hi 2 Eio»(w)i 5 



(2.32) 
(2.33) 

(2.34) 



\n=0 



\n=0 



Moreover one can prove (see [|A|) that if {P n (o;)}^ =0 satisfies ( |2.30|) , then {Q n (cj)}^L 
satisfies it too. Therefore ; 1 < +00. This contradicts (|2. 311). Thus, one has 



(2.35) 



n=0 



meaning that n+ = n_ = 0. We apply Th.VIII.3 from [ |R-S| vol.1. 

ii) If 7^- (0) < +00 then the condition ( 2.35 ) follows immediately from ( 2.31| ). The state- 
ment ii) follows from the statement i). 

iii) One proves it ad absurdum using the previous statement. QED 

Let H re d be some selfadjoint extension of H re d- in the case of type C such extensions are 
parametrized by points of S 1 see [ |R-S| ] whereas in the case of type D the extension is unique. 
Let dE-ji,v (w) be the spectral measure of U re d- 

By the spectral representation of A-n we will call the representation in the space I? (R, dcr-jip) 
where the measure do~n D is given by 



do-R,,V M : = {Q\dEn,v M 0) 



(2.36) 



and the operator H re d acts by multiplication with the identity function on R. Let ?7pjv : 
^red — ► L 2 (R, da-jip) denote the intertwining operator for these two representations. Be- 
cause the polynomials {P n (cj)}^Lq form an orthonormal basis in L 2 (R, da-iip), it is conve- 
nient to write the intertwining operator using the physical notation of Dirac 



U PN = P n (w) ® (n| . 



(2.37) 



ra=0 
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The convergence is understood in the sense of weak topology. The relation of the holomorphic 
representation and the spectral representation is given by the isomorphism of Hilbert spaces 

U PZ : L 2 CD, dfm) L 2 (M, da n ,v) (2.38) 
and Upz one can be written as 

CO n 

fa °£ fiM Vi,)'.«w (239) 

where < , = > is an orthonormal basis in L 2 CD, du-jz). Similarly let 

\ VW( 9 )...W(</») J n=Q V ^ J 

Uzn ■ Hred — > ^ 2 (B, c^w) (2-40) 

be given by 

00 n 

t/ ZjV : = V Z = ® (n| . (2.41) 

^ V* (?)•■■ ft (9") 

We are thus getting the following commutative diagram 

Upz Uzn 

/ \ (2.42) 

L 2 {R,da n ,v) ^ L 2 0(D,^) 

of Hilbert space isomrphisms. 
If the series 

00 n 

^■''^■'" W-.ttOT (2 - 43) 

is pointwise convergent for all a; G [a, 6] and z£D, the isomorphism C7 pz can be represented 
as the integral transform 

(Upz<p)(u)= I V(u,z)<p(z)dim(z,?). (2.44) 



The kernel V (oj, z) of this transform satisfies the 7£-difference equation 

(w - z) V (u, z)=V (Q) V (u, z) + d n V (to, z) . (2.45) 

Of course, from the orthogonal polynomials theory point of view V (uj, z) is nothing else than 
the generating function for the family of orthogonal polynomials under consideration. 

The function V(uj,z) is useful to calculate many important physical quantities of the 
system described by the Hamiltonian H re d- 
First of all, note that using ( 2.44| ) and ( |2.45| ) one obtains 



(v \H? ed z) = [ V(u, v)co n V {u, z) da n p {uj) = ( z + V{Q) + d n ) n Exp n (vz) . (2.46) 
Putting in particular n = we have 



Ex Pn (vz) = / V (u, v)V (u, z) da n ,v (w) (2.47) 
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which is the integral representation of the ^-exponential function Exp-j^ (vz). This function 
satisfies the equation 

d-jiExp-ji (vz) = vExpn (vz) . (2.48) 

Let us recall that Expn (v-) G L 2 (D, dfin) is the expression of the coherent state in holo- 
morphic representation. 

The evolution operator U (t) := e lHredt acts on the function from L 2 (M, do-ji.v) as multi- 
plication by phase factors: 

U (t) tp (w) = e iwt tp (w) . (2.49) 
It enables us to calculate the transition amplitudes between coherent states : 



(v\U(t)z)= / V (uj,v)e iujt V (io,v)da 1ZtV (uj) . (2.50) 
Jr 

The vacuum - vacuum transition amplitude is also important for physicists. It is given by: 

{0\U (t) 0) = I e^da n ,v (w) = T (2.51) 

where fj, n is the n-th moment of the measure dcrup . 

The above shows that the measure da-ji t x> plays a significant role in the description of our 
physical system. The construction dan t D is one of the most important problems which has 
to be solved in order to recover the dynamics. 

In order to give an example see 0] of the solution of this problem, let us recall the notion of 
a simple symmetric operator. 

The Hilbert subspace TLi C Ti is a reducible subspace of the linear operator T : TL — > Ti if 
Til and TL2 '■= Tij^ are invariant subspaces for T and for orthogonal projection IIi : TC — > Ji\ 
one has IIi (Dt) C Dt- The symmetric operator T is simple if there does not exist the 
irreducible subspace of T such that Tj Wl has selfadjoint extension in Hi- 

If the Jacobi matrix of the reduced Hamiltonian H re d is of type C then the series 

00 00 

^|P„H| 2 & ]T|QnM| 2 (2.52) 

n=0 n=0 

are almost uniformly convergent on C. Thus the functions A (iv) , B (uj) , C (to) and D (u) 
defined by 



00 



A(oj) = oj^Q k (Q)Q k (uj), 

k=0 

00 

B(lo) = -l+ujY,Qk(0)Pk(uj), 

k=0 

00 

C(lo) = l + uY,Pk(tyQk{u), (2-53) 

fc=0 

00 

D(lo) = co^2P k (0)P k (io) 

k=0 



are entire functions. If in addition H rec i is simple and closed operator then the spectral 
measure dE-ji,v (w) of its arbitrary selfadjoint extension H re d is localized at the nulls Ui, i = 
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1, 2, . . ., of the function q (u) = B (u) t — D (u) for some tsR. The steps fa, i = 1, 2, . . . of 
the measure da-ji^ (w) := (0| di^x) (w) 0} satisfy the following conditions 



From the identity 



V — ; ^ ? < oo (2.54) 

^^(l + ^l^^)! 2 

oo - 



|n) = P n (tf red ) |0> (2.56) 



it follows that if the Jacobi matrix of H re d is of type D then H re d has simple spectrum. 
Conversely, one can prove (see [^t|) that every selfadjoint operator H = H* with simple 
spectrum may be represented in some orthonormal basis by the formula ( [2.22 ) where the 
Jacobi matrix is of the type D. 



Using (2.3-2~5|) we can associate with it the algebra A-r. with the proper structural function 



1Z. This fact indicates that the algebras of this kind are important tools to investigate the 
symmetry structures of the physical systems with dynamics generated by Hamiltonians with 
simple spectrum. 

In the next section we will describe the situation when neither coherent states nor the 
kernel V (ui, z) do exist. 



3 The integrable systems related to q-Hahn's polynomials 

We will integrate quantum systems related to the algebras A-n and Hamiltonian H re d =: Hab 
with the corresponding structural functions of the form 

K AB (x) = (l-q)x [d gr] (x) - (x) 8 q P (x)] (3.1) 
V AB (x) = (l-q)xd q p(x). (3.2) 



where 



a ( T \ - 1 ( 1 ~ x ) K Q ( 1 ~ g) ~ bl ) x + M Ci-i) 
1 1 " (l- g )[(a 1 (l-g)-6 2 )x 2 + M 21 ' 



f [(a (1 - q) - h) x + b lQ 2 ] [(a (1 - q) - 6Q x + big] 
V [X) 1 [(ax (1 - q) - b 2 ) x 2 + b 2 q 2 ] [( ai (1 - q) - b 2 ) x 2 + b 2 q 3 } + 1 " ) 



+ 



&o(l-g) \ q 3 (l-x)(l-q- 1 x) 



{ ai (l-q)-b 2 )x 2 + b 2 q 3 J (l- q ) 2 (l + q ) 

These functions depend on five real parameters ao,oi and bo,b\,b 2 . It will appear later on 
that it is natural to introduce the following two polynomials: 

A{oj) = aiu + ao (3-5) 
B(u) = b 2 u 2 + biuj + 6 . (3.6) 

It is clear from (37J[37J|) that the pairs of polynomials (A (cj) , B (u>)) of degree one and two 
respectively, taken up to common overall real factor c / 0, parametrize the models under 
consideration. The only condition we will impose on the pair (A (w) , B (w)) is the one given 
by Hab (q n ) > for any n G N. 
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Analogously to the theory of classical orthogonal polynomials there is a g-difference equa- 
tion (an analog of Pearson equation |]Su| ) : 



d q (gB) (uj) = (gA) (uj) (3.7) 

associated with the pair (^4 (uj) , B (a;)). 

We will look for the solutions g (to) of fl3.7| ) which do satisfy the boundary conditions 

q (a) B (a) = g (b) B (b) = 0, (3.8) 

for some fixed a, b such that — oo < a < b < oo. We thus have so called Pearson data 
(A (uj) ,B(u)) on the interval (a, b) C R. 



Proposition 3.1 Let g(uj) be the solution of the q- Pearson equation (3. r i) defined by (A (uj) ,B (a;)) 
and satisfying ( \3. <jj ). Then g( k ' (uj) := g (q k uj) B (quj) . . . B (q k uj) is the solution of Pearson 
q-equation (3.1) associated with the pair 

A^(uj) := q k A{q k uj)+ l ~ i q ^ k d q B{u) (3.9) 

BW(w) := B{uj). (3.10) 
where k € N. If g(ui) satysfies the boundary conditions (3.t) then g^ k \uj) satysfi.es (\3.<fy too. 



Proof : By straight forward calculation. QED 

Let I? ([a,b] , duAB.) be the Hilbert space of square-integrable functions with respect to 
the measure 

do~AB (jjj) = Q (^) dqUi (3-11) 

where 

oo 

d qU j = ^ (1 - q) q k bS (u - q k b) - a5 (uj - q k a?j duj (3.12) 
fc=0 

is the Jackson measure on the interval [a, b]. It will be assumed that the weight function 
g(uj) does satisfy the Pearson g-equation ( |3.7| ) supplemented with boundary condition ( |3.8| ). 

Applying the orthonormalization procedure to the monomials {uj n }^_ C L 2 ([a, b] , doAB) 
we obtain the system of orthonormal polynomials (OPS) {P n (uj)}^ =0 

/ P n (uj) P m (uj) d(JAB M = 5nm (3.13) 
J a 

which is uniquely determined by the pair (A (uj) , B (uj)) and interval. The polynomials 

Pn (w) \ 

J n=0 

the monic OPS associated with {P n (w)}^ =0 . (e.i. P n (uj) = ^-P n (<^) where a n is the coeffi- 
cient of the higest power in P n (w)). 

Theorem 3.1 If\ P n (uj) \ is the monic OPS corresponding to the Pearson data (A (uj) , B (uj)) 

I J n=0 

then the family of polynomials: 

i 



[n] [n — 1] . . . [n — k] 



d«P n (u)\ (3.14) 



ra=0 



where 

l-q k 

forms the monic OPS corresponding to (AS^ (uj) ,B^ k ' (uj)) , with the same boundary condi- 
tions 
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Proof : For k < n — 2 we have 

rb 



P n (uj)uj k A(uj)Q(uj)d q uj = 0. (3.15) 
Using Leibnitz rule, ( ^7| ) and ( |3.8|) we obtain 

= / P n (uj) cu k (d q Bg) (uj) dgu = 

J a 



P n (uj)uj k B(uj)Q(uj) 



b 



d q ( P n (uj) io k ) B (quo) q (qui) d q uj 



d q (P n (uj) (B (uj) - (1 - q) ujA (uj)) g (uj) d q uj = 
b _ 

d q P n (u) (qoj) k (B (uj) - (1 - q) toA (uj)) q (uj) d q uj - 
P n (uj) [k] uj k - 1 (B (uj) ~(l-q) ujA (uj)) q (uj) d q uj. 
The degree of polynomial 

[k}uj k - 1 (B(uj)-(l-q)uA(uj)) (3.16) 
is k + 1 < n, and from ( |3T9 - |3 . 1 0[) one finds that: 

d q P n (uj) uj k g (i > (uj) d q uj = (3.17) 



J a 



for k < n — 2. This shows that the polynomials |j^y<9<j-fn form a monic OPS for the 

Pearson data (A^ (u),B^ (uj)) given by (|3.9H3.10D . QED 

Since, for the rational function IZab (x) of (3.1) one has TZab (0) = the Proposition 
2.1 and it's consequences described in the previous section imply that in the case under 
consideration the following statements are true: 

i) the operator Hab is selfadjoint and has simple spectrum; 

ii) the coherent states do not exist for A-n AB 

hi) the Hilbert space TL re d is unitary isomorphic to L 2 ([a, b] , do~ ab ) ■> with the isomorphism 
given by ( 2.37 ). 

Hence, the Hahn's polynomials {_Pn(^)Yn=Q f° rm an orthonormal basis in L 2 ([a, b] , do~AB)- 
The measure doAB is the expectation value of the spectral measure dEAB in the vacuum 
state |0). 

It is then clear that the properties of Hahn's polynomials are crucial for better under- 
standing of the physical systems corresponding the algebra A-r. ab - The proposition below 
gives a description of some important properties of these polynomials, namely 

Theorem 3.2 (Hahn) Fix a Pearson data (A(u), B(u)) on the interval [a,b] C R. Then 
the following statements are equivalent: 

A. The family of polynomials \p n (uj)X forms the monic OPS with respect to (A(oj) ,B (uj)). 

I J n=0 
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B. The polynomials are given by formula of Rodriques: 

1 



Pn(u>) 



q{uj)B{uj)B ...B [q 



-(n-l) 



(3.18) 



n£N, where c n is a normalization constant. 
C. The polynomials \p n do satisfy the following q-difference equation (Hahn equa- 

I J 71=0 

tion) 

[A ( W ) d g + 5 (w) dgQ^dq) H (u) = \ n K (w) 

u>/zere 



A n = oi[n]+6 2 [n][n-l]<z-( n " 1 J; n = 2,3, 
Ai = oi. 



(3.19) 
(3.20) 



D. Every polynomial of the system I P 1 



}oo 
n=0 



is given by 



n-l 



^ m = n 



,(fc) 



fc=0 ' - 6 2 [— n + l + k] 
where the linear functions 



A<® (u) + B (w) c^Q" 1 



(fe) . (fc) 
u + a 



A^ 1 ) (w) + B (a/) 

(3.21) 

(3.22) 



are defined in (S.i). 



are related by the three term recurrence 



E. The polynomials of the system 
formula 

P n+1 H + K AB (q n ) P n -i M = (w - Pab (q n )) Pn (w) 
wit/j i/ie initial condition Pq (u) = 1. 



(3.23) 



The proofs of the equivalence of A,B,C may be found in the original paper of Hahn 
Q. The recurrence formula ( |3.23| ) is considered there without specification of the form of 
the structural functions TZab^ab- The only assumption made is that they are rational 
functions of the parameter q. A complete proof of this theorem is given in the Appendix B. 

From the considerations of Section 2, it follows that the problem of integration of the 
multiboson system described in 1? ([a, b] , oVab), is reduced to the construction of the measure 
do ab- According to ( 3.11 ) the measure do~AB is given by the density function q(uj) which 
is a solution of the g-difference Pearson equation ( |3.7|) . Let us therefore present all possible 
solutions of ( |3.7| ) from the class of the meromorphic functions. Using (3^-3.6) we can rewrite 
Q3.7|) in the form 



B{qu) 



B(u) - (1 - 



■g(qu) 



b 2 q to + biqu + b 



(6 2 - (1 - q) oi) to 2 + (bi - (1 - q) a ) u + 6, 



■g(qu) 
(3.24) 

and after standard calculations we obtain the following classes of solutions depending on the 
values of the parameters b 2 , b\, bo, a%, ao: 

Proposition 3.2 One has the following subcases of the solutions of the q-difference Pearson 
equation (3.1): 
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i) and 6 2 — (1 — q)a\ / 0, then 



q{u) 



{ q a i g)oo( 9 fr i g)oo 
00 



(3.25) 



where 0/ 0, 6/0 are roote 0/ t/ie polynomial B(u) and c / 0, d / are roots 0/ t/ie 
polynomial B(uj) — (1 — g)uM(cj). 



ii) If bo and b\ — (1 — g)ao / and 62 — (1 — g)ai = , taen 

(lT> 9)oo(^"; 9)00 



(3.26) 



where a / 0, 6/0 are roots 0/ t/ie polynomial B(ui) and c / is t/ie root 0/ the 
polynomial B(oj) — (1 — q)ujA{uj). 



iii) 7/ 60 / and 61 — (1 — g)ao = and 62 — (1 — g)ai = , then 

qu qu 

Q{u) = (— ;g)oo(-^-;o)oo, 



(3.27) 



where a / 0, 6/0 are roots 0/ t/ie polynomial B{uj). 
iv) 7/ 60 = 0, 61/O and 61 — (1 — g) ao / and 62 — (1 — g) ai / and 62 / , t/ien 



o(w) = uf' 



y r i 9)oO 



(3.28) 



where a / is t/ie root 0/ t/ie polynomial B{oj), c / is t/ie root 0/ t/ie polynomial 
B{uj) — (1 — g)ctM(cu) and g~ r 



6i-(l-g)a 



v) 7/ 60 = 0, 61/O and 61 — (1 — g) ao / and 62 — (1 — g) a\ = and 62 / , then 

(3.29) 



where a / is the root of the polynomial B(uj) and q 



qbi 



6i-(l— q)ao 



vi) If 60 = 61 — (1 — g) ao = 0, 61 / 0, 62 / and 62 — (1 — g)ai / 0, then 
a) 

q(lo) = uj t 



(f;< 



(-u;q) 00 (-qu> 1 ;g) c 



for q~ 



>0; 



(3.30) 



b) 



g(co) = uj r 



( 



qu) . 



(w;g)oo(gw x ;g) c 



where a / is the root of the polynomial B{oj) 



(3.31) 
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vii) If bo = bi = 0, b\ — (1 — q)a§ 7^ and 62 / 0, then 
a) 



/ \ r (-a;;g)oo(-ga; 1 \q) 00 . . 

= w 773 — r (3.32) 



for <T = <^^r > 0; 



b) 



= u — — (3.33) 

I c 1 Q)oo 



where c 7^ is the root of the polynomial B{oS) — (1 — q)u)A{uS). 
viii) If 60 = bi = b\ — (1 — q) a$ = 0, and 62 — (1 — 9)01 7^ and 62 7^ 0, t/ien 

q{uj) = oo r (3.34) 



/or q- 



62— (l-g)ai 



Proof : The subcases i), ii) and iii) are easily obtained by iteration. The points iv)- viii) 
are proved by calculation of Laurent expansion coefficient and application of the Ramanujan's 



identities (see |G-R|1 ). QED 



We can now determine the interval of integration in ( 3.1 3| ) and determine the conditions 



on polynomials A{u) and B{uj) such that the measure do~AB is positive (ie R{q n ) > for 
n 6 N). It will be convenient to express the conditions on A(u) and B(u>) in terms of roots 
of the polynomials B(u) and B(u) — (1 — q)ujA{u). 



Proposition 3.3 The measure do~AB is positive and the condition (3.8) is fulfilled if and 
only if (in the notation and classification of Proposition 3.2) 

i) The integration interval is [a, b] with a < < 6 and c, d satisfies one of the following 
conditions: 

a) c = d, 

(3) c < a and d > b, 
7) c, d < a, 

5) there exists K £ N such that q K ~ l a < c,d < q K a, 
e) c,d> b, 

C) there exists K £ N such that q K b < c,d < q K ~ l b. 

ii) The integration interval is [a, b] with a < < b and c < a or c > b. 
iii) The integration interval is [a, b] with a < < b. 
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iv) This case splits into two subcases: 

1) For a > the integration interval is [0, a] and c < or c > a. 

2) For a < the integration interval is [a, 0] and c < a or c > and r has to be 
such that a r > 0. 

v) This case splits into two subcases: 

1) For a > the integration interval is [0, a]. 

2) For a < £/ie integration interval is [a, 0] and r aas to 6e such that a r > 0. 

vi) This case splits into two subcases: 

1) For a > the integration interval is [0, a]. 

2) For a < i/ie integration interval is [a, 0] and r aaue to 6e such that a r > 0. 

vii) In this case IZAB(Q n ) ar & not positive for n large anough. 

viii) In this case 7ZAB(Q n ) = and T>AB(Q n ) = for n G N. 
Proof : 

i) The equation B(uj)q(uj) = is solved by aq~ k+1 and by bq~ k+1 for £ H. For any 
function f{u) and any GN, using (ETIID and ( [Up one can obtain 

a g-fc+l pbq~ k+1 

f{u)duAB{u) = =/ f(oj)da A B(u) (3.35) 

aq-i+i Jbq- l + 1 

and 

f{uj)d<JAB{u) = / f{<jj)d(jAB{u). (3.36) 
Hence have the integration interval is [a, 6]. The condition of positivity of d(7AB{^) 

f f{u)doA B {u) > /or/ > (3.37) 

is equivalent to 

aa(<fa) < and bg{q i b) > /or i = 0, 1, (3.38) 

The continuity of g at uj = gives a < < 6 and the inequalities 

a(g i a)>0, g{q i b) > /or i = 0, 1, . . . . (3.39) 

which are solved by o;)-C)- 

The proofs of ii)-viii) are similar to the one above. QED 

The above class of orthogonal polynomials, which we call the q-Hahn polynomials, con- 
tains, as a special cases the families of orthogonal polynomials well known from literature. 
Using a very good paper |K-S| we obtain the following identification: 

1) Putting in \)j3) d = 1 we have the Big q-Jacobi polynomials. If in addition we put 
b = q and c = | we obtain the Big q-Legendre polynomials. 

2) Putting in ii) b < 1 and c = 1 we obtain the Big q-Laguerre polynomials. 

3) Putting in iii) b = 1 we obtain the Al-Salam-Carlitz I polynomials. If in addition we 
assume a = — 1 we obtain the Discrete q-Hermite I polynomials. 
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4) Putting in iv)l) a = 1 we obtain the Little q-Jacobi polynomials. If in addition we put 
c = - and r = we obtain the Little q-Legendre polynomials. 

5) Putting in v)l) a = 1 we obtain the Little q-Laguerre/Wall polynomials. 

6) Putting in vi)l) r = we obtain the Alternative q-Charlier polynomials. 
We will find now the equations for the moments 

Hn= f u n da AB (uj) (3.40) 

J a 

of the measure do ab- 

From Section 2 it is clear that once the moments are known one may determine many 
important physical characteristic of the system under consideration. 

Multiplying (/-difference Pearson equation ( |3.7|) by u n q n and using Lebnitz rule for q- 
derivative we obtain the following three-term recurrence equation 

- N (hfhi+l + h^n + &oMn-l) = q n (oi/in+i + a fj, n ) (3.41) 

for n > 1, and 

aini + a Q no = 0. (3.42) 

The initial rule /io = X? doAB^) for this recurrence can be calculated in straightforward way. 
In terms of the notation and classification introduced in Proposition 3.2 we have 



i) 



(?; ?)oo(g|; q)°o{ql;q)oo{^q) 



Mo = (1 - q)(b - a) (3.43) 

(|;ff)oo(S;ff)oo(J;ff)oo(i;?)c 



ii) 



n \ / j , (q;q)oo(qa;q)oc(qr><i)°o ,o A a\ 

Mo = (l-g)(6-o) ° fc (3.44) 

(f,q)oc&q)c 



I oc 



111 



IV 



v) 



/i = (l - q)(b- a)(g;9)oo(g-;g)oo(gr;9)oo (3.45) 

a o 



Mo = (l-g)a r+1 7^f- (3.46) 



Mo = (l-g)a r+1 (g;g) r (3.47) 
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vi) a) 



-i(q;q)oo{-aq r+1 ;q) 

Mo = K 1 ~ <l)u 

b) 



Mo = (1 - q)a T -j; . g r 3.48 



, r +i (g;g)oo( Q g r+1 ; g) 



/io = (l- 9 )a r+1 ^^' (3-49) 



oo 



Let us note that replacing in iv)-vi) r by r + n, n £ N we obtain the moments for the 
corresponding cases. 

In order to consider the cases i)-iii) let us introduce a real function fi satisfying the 
equation 

(l-w)B (Q) m (w) + (1 - q) ojQA (Q) n (u) = 0. (3.50) 

It is easy to check that fi (q n+l ) satisfies the recurrence equation (|3.41| ), hence ^{q n+1 ) = Mn- 
Reexpressing ( |3.50| ) in the form 

B{ - Q) u)fx (w) = (3.51) 



B(q-iQ)-(l-q)q~lQA(q-iQ) 
and the equation ( |3.7| ) in the form 

B(qcu) 



Q]q(u)=0 (3.52) 



B(uj)-(l-q) loA (w) 

one may observe some symmetry between the equation on g and the equation for the moment 
function jj,. After the substitution of the form Q — > qui and u — > Q the operator from ( |3 . 5 1] ) 
transforms into the one of ( |3.52| ). The equation ( 3.51 ) as well as the equation ( |3.52|) can be 
easily solved. 

For example, if we assume that B(l) = then ( |3.51 ) can be written in the form 

&r,h{u) = n(u), (3.53) 
where &n is ^-derivative. The function 1Z is here given by 

, x B(u) 

v ' B(q~ l uj) - (1 - q)q- 1 ioA(q- l uj) v ' 

Then one of the two linearly independent solutions of ( |3.51| ) is simply ^-exponential Exp-ji- 
In this case it is given as the basic hypergeometric series 

Mi (uj) = Exp n {u) = 3 $2 ( a i Vl (3.55) 



where a / 1 is the root of the polynomial B(u) and c, d are roots of the polynomial 
B(q~ 1 uj) - (1 - q)q~ l ujA{q- l uj). The function 3 $ 2 is defined in gj§ , |^Sj . 
The second solution ^2 (k>) is related to Ml (w) by the following formula ((/-version of Wron- 
skian) : 



„A(«9;g) c 



M2 M Mi - M2 (<?w) Ml (w) = x - — ^-r-^ (3.56) 

{^q)oo 



where q x = ^, a - H^hl^ 



bo ' hi 
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Any solution \i (to) is a linear combination of \i\ (cj) and /x 2 (uj). We are then getting the 
following formulae for moments 

/in = li {q n+1 ) = cx/xi (g n+1 ) + c 2/ x 2 (<f + 1 ) (3.57) 

where the constants ci,c 2 are determined by 



Mo 



/ do-^B = c x /xi (g) + c 2/ u 2 (g) (3.58) 

J a 

ao (ciMl (?) + C2M2 (<?)) + «l (ciMi (<? 2 ) + c 2 /i 2 (g 2 )) = 0. (3.59) 

Appendix A 

The affine difference calculus and q-Hahn's orthogonal polyno- 
mials 

In this section we present the preliminary considerations related to the calculus generated by 
the action of the affine group A + on the real line. Let us define the linear representation of 

(£ q , h <p)(x):=ip(qx + h), {q, h) G A + = {(q, h) : q > 0, liGK) (A.l) 

acting on the functions 99 from the algebra 5- Since our consideration will be formal in its 
character we do not impose any additional conditions on 
According to |h1 we introduce the derivative operator 



<p(x) — (f (qx + h) 
x — {qx + h) 



(*wo (*) : = y - "i- r.T ; ( A - 2 ) 



as a natural generalization of the q- derivative d q := d Qj o and of the difference derivative 
dh ■= di,h- 

The Leibnitz rule for the derivative d q ^ is 

(dq ;h cpifj) (x) = (d q>h <p) (x) ip (x) + (£q th <p) (x) (d qjh il>) (x) (A.3) 
There is also the following equivariance property: 

£c,t ° d Q,h O £ c ,t = Cdq,ch+(l-q)t (A.4) 

and enables us to reduce (q, h) —analysis to q— analysis. We have for example 

dq, h = £; l h o9 ? o£ iJl . (A.5) 

L >l-q '1-9 

Let us now solve the equation 

dq,h<P = Q ( A -6) 
for the given function g G In the order to do this, we apply the operator Z k h to 



and we find that 

ifc+l 



^ (x) - £*+V (x) = q k [(1 -q)x-h) £* >h Q (x) . (A.7) 
Summing up both sides of the identity dA.7|) with respect to k we get 

if (x) - ip (x°°) =Y,[x-(qx + h)} q k Q ( q k x + ~^~h j (A.8) 
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where 

x°° = lim ( q k x + ^—^-h] = . (A.9) 

fc^oo y 1 — q J 1 — q 

The equation (|A.8|) , justifies the following definition of the (q, h) —integral 

f Q(x) = T Q (t) d q , h t := V [x - (qx + h)\ q k Q (q k x + ^^h) . (A.10) 

Jq,h Jx<™ ^0 ^ ' 

The (q, h) —integral operator is the right inverse of the (q, h) —derivative operator 



d q , h o / =id, (A.ll) 

Jq,h 

and, moreover, 

/ °d qt h = id- <5i,oo- ( A - 12 ) 

J q,h 

The oprerator 5\ i00 is an idempotent operator defined by 

(5 1 , 00 <p)(x) = <p(x co ), (A.13) 

projecting the function on the constants. 
Like in ( |A.E| ) we have 



= ST 1 j! _ o o^j^, (A.14) 

q,h '1-8 

which reduces (by the translation authomorfism £, the (q, h) —integral to the Jackson 

'9-1 

integral J q = j q0 - The integration on the interval [a, b] can be defined by 

pb pb pa 

/ Q{t)dg !h t= g(t)d q>h t- g{t)d qA t. (A.15) 

If q — > 1, the calculus presented above corresponds to the difference calculus. For h — > 
one obtains g— difference calculus. The differential calculus will be obtained when q — > 1, 
/i -> 0. 

Let us mention finally that the identities (|A.5| ) and ( |A.14 ) enable us to reduce (q,h)- 
calculations to the g-calculations. This property motivates us to discuss the case of the q— 
analysis in this paper. 



Appendix B 

Proof of the Theorem 3.2 

The monic OPS {P n (w)}^L ^ s uniquely defined by the weight function q{uj) and the 
interval [a, b]. In order to prove the equivalence of the properties A and B it is sufficient to 
show that the system of polynomials defined by ( 3.1§| ) is a monic OPS. In order to do that 
let us reexpress the function 



q(uj)B(uj)B (q^uj) ...B 



q-^oo 



(B.l) 



in the following way 

F k (lo; n) = q(lu)B (to) B (q^cu) ...B (qr^ -1- *^) R Kn (to) k = 0, 1, . . . , n - 1 (B.2) 
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F n (u; n) = g {lo) R n , n (u) 



(B.3) 



where Rk, n {lo) is a polynomial of degree not greater than k < n. These polynomials do and 
satisfy the recurrence formula 

R k+1 , n {lo) = A {lo) R k , n {qu)+B (q^^u) d q R k , n {u)+—^ - —R k ,n {qui) 



(l-q)u 



for k = 0, 1, . . . , n — 1 with the initial condition Rq )TI {lo) = 1. 
For k < n, applying (|B.2j) we have: 



(B.4) 



fc(fc+i) , ~ /"" fc(fc+i) , 

2 u*P n (u) g(u)dgU = c n q 2 uj K d™F (uj;n)d. 



fc(fc+l) , 

c„a 2 uj (u)\n 



c n [k] 



fc(fc+i) 



W^ufQ {lo) B {lo) B (q^u) ...B (q-^ 1 -^- 1 )^) R k _ 1>n (uj) " - (B.5) 



-c n [k] I q^u^d^Fo^^d^ 



(_!)* [fc] [ fc - 1] . . . [1] Cn / d n ~ k F {lo; n) d q oo 



{-l) k [1] . . . [k] g {lo) B {lo) B {q- l u) ...B (q~^- 



(fe-1)) 



lO 



0. 



This shows, that polynomials P n {lo) , n E NU {0}, form an OPS. By the proper choice of the 
normalizing constants c n one can obtain the monic OPS. 
B => C 

We have proved the validity of the q— Rodrigues formula for any Pearson data. So by 
Theorem 3.1, we have 



1 



— d q P n {lo) = cji 



1 



-0, 



,n— 1 



? « {lo) ( w ) ^(1) . ..£?(!) (g-<"- 2 >w)] (B.6) 

for the Pearson data (A^ {lo) , B^ {lo)) given by ( |Q[jriC| ). Using now the equality 

Qg {lo) B {lo) ...B (V (n ~ 1} ^) = M (w) (g- 1 ^) (V (n ~ 2) ^) (B.7) 

and substituting (B.6) into ( p,18| ) we find 



Pn {lo) = <*JTj& [Q- 1 ^ M («) B W (^ (n " 2) ^) 



(B.8) 



. 9 -(»-l)l Q-l#»-l L« ( w)jB (l) (cu)jB (l) fg-("-2)^ 

o(a0 y L v \ / 



c n q 



g{uj) 



-n-1 



<r (n ~ 1} [n] U (w) <9 g P„ (w) + 5 {lo) 8 q Q- l d q P n {lo) 
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We have proved ( |3.19| ). In order to show (3.20) we compare the coefficients of x n on both 
sides of ( 3.19| ). Additionally we find the formula 



-n-l 



[n\ q 



(B.9) 



for the normalizing coefficients. Later we will use ( p.9[ ) for the calculation of c n . 
C^B 

The proof goes by the induction. It is easy to see that for n = 1 ( p. 18 ) follows from ( |3. 19| ) ■ 
Let us assume that it is true for n — 1. We prove it for n: 
From our assumption and Theorem 3.1 we have 



n 



-d q P n (uj) 



i^TT^ 1 (e W H B {»)... B (q-^u 



(B.10) 



Using ( 3.18j ) and the Proposition 3.1, we obtain thesis after a simple calculations. 
C^D 

From (ET2CD we have 



P n ( W ) = M (A ( W ) + B M fyQ" 1 ) i_0 ff P n (, 



n 



(B.ll) 



According to the Theorem 3.1, the polynomials | ^d g P n form the monic OPS with 

respect to the data (A^ (uj) , B^ (uj)) given by (pH3~T(i|). We can thus apply the formula 
( B.ll| ) with (A^ (uj) , B^ (u)) to the polynomial -^d q P n (uj). Repeating this procedure 
n-times and using the formula 



AW = a[ k) [n] + b 2 [n] [n - 1] q 



(n-l) 



(B.12) 



where is defined by flOSj) , we obtain ( ggg ). 
B => £ 

In order to prove that the recurrence formula ( |3.23 ) holds we use the identity 



d: 



q(uj)B(uj)B {q~ l UJ) —B (q 



-(n-l) 



= q(uj)B (uj) B (q-'uj) ...B (q~^-^ R k>n (uj) 

i I5.13i 

where the polynomial Rk, n ( w ) satisfies the recurrence equation (B.4). From the eq. ( B.13| ) 
and Rodrigues formula one has 



[UJ 



Pn ( w ) — c nPn 

Let us denote by a k , /3j~ and 7^ the three highest coefficients of the polynomial 



(B.14) 



Rk,n (u) = akUJ* + /3 k uj k 1 + 7 fc w 



k-2 



+ 



(B.15) 



After substituting (|B.15| ) into ( |B.14|) and comparing the coefficients of the monomials 
ui k and uj k ~ l we obtain the following system of the recurrence equations 



UJ 



k+l 



(01-62 \-2n + 2 + k])q k a k 



«fc+i : 

P k+1 = (a -b 1 [-n+l])q k a k + (a 1 -b 2 [-2n + 3 + k])q k ~ i p k (BAG) 
7fc+i = b [k]a k + (a -b l [-n + 2})q k - l l3 k + (a 1 -b 2 [-2n + A + k})q k - 2 lk . 
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One can solve them by iteration: 

k 



fc(fc-i) 

1=0 



a k = q^^~[[-l(a 1 - b 2 [-2n + 2 + I]) (B.17) 



fJ k = Mr - a °" 6l r [ " n + 1] 1 -« fc (B.18) 
H q n ~ l a 1 -b 2 [-2n + 2] V ; 

(1 _ q n) (1 _ g n-l) (a() _ 6i [_ n + 2]) (fflo _ &i [_ n + x]) + b(j _ ^ [_ 2n + 2]) 

(B.19) 



7fe (1 - g) 2 (1 + g) ' (a 1 -6 2 [-2n + 2])(a 1 -6 2 [-2n + 3]) 



Thus for the monic polynomial 

P n (a,) = w n + /? (<f ) Z" 1 + 7 (g n ) a;™" 2 + ... (B.20) 
we find that the coefficients 

0(q n ) ■= ^ 
7(9") := J 

are given by the rational functions (|3.3| - [3.4| ), Using the three terms recurrence relation ( 3.24 ) 
and ( |B.15| ) we obtain the formula (3.1-3^) for the structural functions 9\ and D . 

The recurrence formula ( 3.23| ) rewriten for the orthonormal polynomials {P n (w)}^L takes 
the form (|2.2§| ) which means that the Hamiltonian H re( i ( p.22[) has the Jacobi matrix of the 
type D. Thus (see JA[ ) this Hamiltonian is essentially selfadjoint and has simple spectrum. 
This shows that there is a unique measure do ab such that 



\ P n (u) P m (lo) do AB (a;) = 5 nm . (B.21) 

J a 
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